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The problem of two electrons in a two-dimensional random potential is addressed nu- 
merically. Specifically, the role of the Coulomb interaction between electrons on localization 
is investigated by writing the Hamiltonian on a localized basis and diagonalizing it exactly. 
The result of that procedure is discussed in terms of level statistics, and the expectation 
value of the electron-electron separation. 
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Introduction. The effect of the Coulomb electron-electron interaction (EEI) on 
electronic localization in a random potential is investigated by computer simulation 
for a system with two electrons, using level spacing statistics. In the absence of inter- 
actions, it has been shown (see e.g. m) that the distribution shifts from Poisson to 
Wigner as the system goes from being strongly localized to delocalized. For interact- 
ing systems, such correspondence has never been proven. Still, level spacing statistics 
has been used often to assess localization. Here, we attempt to make a connection 
between localization and level spacing statistics for interacting systems. 

The combined importance of disorder, interaction, and elastic tunneling poses a 
very difficult problem. It was suggested twenty years ago that the EEI can delocalize 
the electrons [|2|, but as yet a firm answer is lacking. Computationally, the main 
difficulty is the huge phase space for systems of reasonable size [[| . Existing work 
fC) , (|| - |l(|) resorted to various approximations. The 2-electron problem for reasonably 

large systems can be solved without such approximations, double occupation of sites 
can be accounted for, spin and exchange included, and the entire phase space can 
be examined. The motivation for the problem considered here is the ability to study 
the legitimacy of the approximations made in the finite-density works, and thus shed 
C^ ' some light on this more complex problem. We hope some insight can be gained into 

the mechanisms at play in the experimentally observed metal-insulator transition in 
2D Jll[ . Previous relevant studies on the 2D random system with interactions include 
[|J, finite-size scaling of three and four M, and two fi [lOf spinless electrons. All these 
works have concluded that the interaction enhances derealization. In S a crossover 
from Poisson to Wigner was found, while 0, M reported a sharp transition. 

Computations. The procedure used here differs from methods used in || || [10| 
in that it is virtually exact in a tight binding scheme: spin is not ignored, fluctuations 
in the overlap integrals are not neglected - nearest-neighbor (n-n), next n-n, and 
next-to-next n-n coherent tunneling processes are accounted for. The method we 
use is the following. First, singlet and triplet configurations are written for every 
pair of sites in the system. These configurations are constructed from one-electron 
s-orbitals ip of radius as centered on the sites. The Hamiltonian is written in this 

1 Fax: +1 661 664 2040. email: jtalamantes@csub.edu 

2 permanent address: Institute of Physics, Technical University of Budapest, Hungary 



On 



o 
o 



X 



representation and diagonalized. The resulting eigenstates are analyzed, and the 
level spacing distribution is computed for several system sizes. The basis set used 
are the two-electron configurations 4> a b for a given pair of sites a and b. The 4> a b 
are constructed by symmetrizing or antisymmetrizing products such as (f a {V)ipb{2) of 
one-electron orbitals for electrons 1 and 2. We consider the Hamiltonian 

2 2 

H = Y(T a + V a +e a ) + —, (1) 

^-i Kri 2 

where a labels the electrons, T, V, and e are the operators for the kinetic energy, the 
core potentials, and the random potentials, respectively; e is the electronic charge, k 
the dielectric constant, and ri2 the distance between the electrons. e a is chosen from 
a box distribution — W/2 < e < W/2, with W equal to the n-n Coulomb energy. 

The integrals corresponding to the matrix elements are performed numerically. 
An L x L lattice is set up, and H is diagonalized for the parameters L and r s (the 
n-n distance in units of as). The resulting eigenstates J, and eigenenergies Ei are 
investigated in two ways: (i) the distribution p(s) of nearest-neighbor level spacings is 
obtained. (We dropped 100 states or so from the band edges.) As in g, a parameter 
r\ = (var(s) — 0.273)/(1.0 — 0.273) is computed as a measure of how close p(s) is 
to a Poisson (r/ = 1) or Wigner (77 = 0) distribution; (ii) the expectation value 
A/ = J2ab r ab\Ai,ab\ 2 of the e-e separation is computed from the eigenstates of (ffl). 
r a b is the distance (in units of as) between a and 6, and Ai^ are the coefficients 
defined by I = 2\labAi,ab<t > ab- One expects that in the localized regime A/ is strongly 
correlated with Ej (larger A/ correspond to smaller Ej), whereas A/ should become 
essentially independent of £7 as configuration mixing increases. 

For definiteness, as and k are taken here to be 10A and 3 respectively - these 
values seem appropriate for 2D systems. This choice of parameters yields an effective 
mass m* = 0.16m, with m the electron mass. The charge on the sites is taken to be 
|e|. Cyclic boundary conditions are used. 

Results and discussion. Runs were performed for L = 4, 6, 7, 10 and r s = 

5, 6, 7, 8, 9, 10, 11, 12. In every case the number of "samples" was sufficient to obtain 
~ 1.5 x 10 4 levels. Fig. 1 presents n(r s ) for each value of L; furthermore, to show 
that the clear drop in 77 with decreasing r s is due to the EEI (and not simply weak 
localization due to system size), fig. 1 also shows r\ without the last term in eq. (QJ). 
It is noted that there is no clear small-size scaling behavior, suggesting a crossover 
(r s ~ 9-11), rather than a transition. This is in agreement with H, but differs from 
O 0- T ne differences may be due to differences in models and choice of parameters. 
Fig. 2 presents fits for X(E) with L = 7. It represents a single realization of 
the disorder, but it is typical of other realizations. The numbers on the graph refer 
to the values of r s . The solid lines represent singlets and the crosses triplets. The 
lines towards the bottom right refer to the upper band, which results from on-site 
repulsion. Features worth noticing are (i) As expected, spin plays no role deep in 
the localized regime (since the exchange energy is proportional to (ip a \ipi 1 )' 2 ). Perhaps 
more unexpectedly, spin also seems to have little importance near the crossover to 
delocalization. (ii) Configurations with doubly-occupied sites play little or no role 
near the crossover - they become important only for r s ~ 2. (Hi) In the localized 



extreme (i.e. r s = 12), the rather short span in energy of the eigenstates comes from 
the random and Coulomb energies; as r s decreases, the broadening and shift of the 
range of eigenergies is attributable to the growing off-diagonal energies of H . 

Examining X(E) we observe A to decrease sharply with increasing E for large r s . 
This is easily understandable as an increase in the repulsion energy with decreasing 
A. The dependence of A on E weakens as r s decreases and configuration-mixing takes 
place - A becomes nearly independent of E for r s = 5. The crossover takes place 
around r s = 5-7, which is somewhat lower than the crossover in fig. 1. The value of 
A (for r s = 5) is reasonably close to L/2, the maximum possible separation (except 
in corners) when cyclic boundary conditions are used. This provides evidence of 
collective derealization Q . While the large off-diagonal energy shows the importance 
of elastic hopping, the persistently large e-e separation shows that the electrons move 
in a correlated fashion to stay apart. Studies of X(E) in non- interacting systems show 
a lack of functional dependence, with A varying at random over a large range of values. 
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Figure 1. r\ vs. r 3 for L — 4,6,7,10. For comparison, two plots (L = 6,7) have been 
included for which the EEI is absent. 

Figure 2. A vs. Ei for L = 7 and a given realization of disorder. Solid lines are fits for 
mixtures of singlets (lower bands correspond to doubly-occupied sites), and crosses are fits 
for mixtures of triplets. The numbers indicate the corresponding values of r s . 

Conclusions. The model used here yields zero density in the thermodynamic 
limit, and so no definite claims or comparisons with experiments can be made; never- 
theless, it is interesting that even this model gives electronic derealization for densi- 
ties ~ 10 12 cm" 2 , which is only an order of magnitude off the critical density observed 
(10 11 cm" 2 ) in Si-MOSFETs Q. Possibly a different choice of W might give a better 
value for the critical density; nevertheless, many-electron effects probably also play 
an important role in derealization at finite electron densities. 

Where collective hopping of the two electrons is coherent (r s < 6), A can be 



interpreted as a coherence length. It is of interest to note that a crossover of p(s) 
from Poisson to Wigner occurs at a somewhat larger r s than the crossover from a large 
variation in \(E) to X(E) ~ const. This suggests that (unlike in the non-interacting 
case) for interacting systems, delocalization in real space requires a somewhat larger 
overlap for n-n sites than does the transition to Wigner statistics. We do not yet 
understand this result well, but note that it is in keeping with previous work Jl3[ . 
There it was observed that the wavefunctions are "swiss cheese-like" without EEI, but 
space-filling with the interactions. Thus, while the EEI may make the wavefunction 
extend over more sites, it does not similarly increase its spatial extent; therefore A 
might require a larger overlap than p(s) (for a crossover to take place) because the 
EEI makes the wavefunction more compact in real space. 
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